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Abstract. We analyze the time evolution of an initial spatial coherence for a two level atom whose internal 
degrees of freedom interact with a single mode of a cavity field. When the qubit-field subsystem is taken 
as an environment, the translational dynamics experiences a decoherence process which may be encoded 
in a decoherence factor D. We find that the field statistics affects D through the alternative paths the 
system-environment may follow along their entanglement, while eventual field phase properties give rise to 
an imaginary part of D which is related to the atomic translation. From the decoherence perspective, we 
analyze the relation between the atomic momentum and the imaginary part of the atomic spatial density 
matrix, and some considerations on its asymptotic behavior are brought into question at the conclusion of 
the paper. 

PACS. 03.65.Yz Decoherence; open systems; quantum statistical methods - 37.10.Vz Mechanical effects 
of light on atoms, molecules, and ions 

1 Introduction vation involves ignorance of a subsystem (a part of the 

universe) and this procedure defines the "facts" that can 
Decoherence program [HHSHIMS] may be considered a ^^^^-^^^ ^ quantum system. Accordingly, Landsman 
serious attempt to overcome the basic dilemma of Quan- ^ ^^^^^^^ ^^^^ ^ measurement, a fact or event in QM 
tum Mechanics (QM): Why the world appears as classical, -^^^-^^ non-observation, or irrelevance, of a certain 
despite of its underlying quantum nature that allows for ^^^^ ^^^^^^ question, 
arbitrary superpositions of states? An implication inher- 
ent to this program is the possibility of dividing the world Remarkable experiments on the decoherence of meso- 
into subsystems. As suggested by Zeh [7l|4], any obser- scopic coherent fields in a cat state and complementar- 
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ity experiments on Rydberg atoms interacting with mi- 
crowave cavities have been performed by Brune et al. [8] 
(see also the review of Raimond, Brune, and Haroche [9j). 
Similar foundational aspects of complementarity, which- 
way information and quantum erasure have been studied 
by Storey, CoUett, Walls [11111], Scully, Englert, Walther 
[12] . Storey, Tan, CoUett, Walls [13], Diirr, Nonn, Rempe 
|14j , in systems implying interactions and correlations amonj 
atomic internal and external dynamics and laser or maser 
cavities. 

It is usually understood that the irrelevant subsystem 
works as a reservoir, that is, it consists of the innumer- 
able degrees of freedom of the environment with which 
the system of interest interacts. However, in a natural 
way one may extend the analysis to environments with 
a few degrees of freedom (see for example [T51IT61IT7] ) and 
to pay attention to the role that the density of the states 
the irrelevant subsystem can effectively accede to, plays 
in the decoherence process. The study of these cases in 
which the system may follow a relatively small number 
of alternatives, can be useful for a grasp in the quantum- 
classical transition. For instance, non local correlations of 
entangled Bell's states may be suddenly and quite irre- 
versibly destroyed by an environmental action of a sub- 
system of continuous one-degree of freedom [17j . Coher- 
ences and non local correlations are in fact very sensitive 
to any environmental action, and the ensuing information 
leakage towards the irrelevant part of the system is at the 
origin of the classical landing. 

Decoherence effects on qubits caused by the atomic mo- 
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tion and, viceversa, effects of the qubit-field interaction 
on the atomic spatial coherences have been considered in 
different contests (see for example [HdTHHITllgnilgllllll 
[231124]). In particular, in Ref. [23]l24] the optical Stern- 
Gerlach (SG) model was used to analyze complementarity, 
which-path information and quantum erasure. Using the 
same model, the time behavior of the entanglement be- 
;tween the internal dynamics of a two- level atom, the zero 
point cavity field and the transverse translational variables 
of the atomic center of mass, was analytically studied in 
Ref. [10] ■ Disregarding the translational dynamics by trac- 
ing on the relative variables, the field-qubit coherences go 
to zero at long time. The diagonal form asymptotically 
attained by the reduced density matrix, clearly indicates 
how the field-qubit system goes towards the separability. 
In addition, if one choose the qubit as the system of inter- 
est, that is, tracing also on the field variables, one trivially 
recover the classical diagonal form of the qubit density ma- 
trix [25] . These are two examples of decoherence in which 
a decisive role is played by the irrelevance of a single de- 
gree of freedom, with a continuum of accessible states, 
that is, the conjugate variables that describe the atomic 
kinematics along the cavity axis. 

In the present paper we will inquire on the atomic 
spatial coherences caused by the qubit-field interaction. 
Using the same model, we analyze the time evolution of 
an initial spacial coherence under the effect of the in- 
teraction between the field of an ideal cavity and the 
atomic internal dynamics. Looking at the qubit-field vari- 
ables as the irrelevant part of our system, a decoherence 
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effect of the atomic position follows, and we show that consider Rydberg atoms interacting with microwave cav- 
the density matrix that describes the translational dy- ity field. The same kind of interaction is used by Scully 
namics at a generic time t may be simply factorized in et al. 12J in their analysis of complementarity, based on 
terms of the initial density matrix and a decoherence fac- matter- wave interferometry. The possibility of realizing 
tor, p(x, x'; t) = p{x, x' \ 0)D{x, x'; t), where D{x, x'] t) will initial coherent spatial distributions using Rydberg atoms 
depend on the field statistics. We analyze the decoherence in microwave cavities, has actually been called in ques- 
process for different initial configurations of the qubit-field tion in Ref. On the other hand, strong-coupling con- 
subsystem. Starting from the more general state (pure or ditions readily obtained for these atoms and high qual- 
mixed) of the cavity single mode, we then specialize to ity factors Q, allow to neglect, with some accuracy, both 
some cases, including the thermal and coherent states, spontaneous atomic emission (tq ^ 10~^sec) and cavity 
We also give the results for the so-called eigenstate of the loss (tc ^ 10~^sec) [ffj for the atomic flight times we will 
Sussking-Glogower phase operator p51[?fll^ . In all these use to obtain full decoherence (Tq < 10^^ sec). 
cases we find a loss of coherences of the atomic position, 

similar to the well known phenomenon of the Rabi oscil- 2.1 Model and initial configuration 

lations collapse in the usual Jaynes-Cummings model [291 

The optical SG model is particularly suitable for an ana- 

[30l[3Tll32l[33l . Because of the discreteness of the variables 



lytical study of our subject. It consists, as known, of the 

that account for our environment, a partial revival of the 

usual Jaynes-Cummings model in which, however, the dy- 

coherences follows at relatively long times. 

namics of the atomic center of mass along the cavity axis 

Finally, we will address to the question of the decoher- 

is taken into account, 

ence process in the presence of a mean value of the atomic 

momentum. We find that for an increasing mean atomic H = \-huj ( a^a + Sz + - ]+hekx (a^ S-+a S+), (1) 

Zm \ I ) 

momentum, as it is for some configurations here consid- 

and correlates with the dynamics of the other subsystems. 

ered, the imaginary part of p(x, x' \ t) does survive in some 

This Hamiltonian describes, in the rotating wave approx- 

regions of the plane {x,x'). 

imation, the resonant interaction of a two-level atom of 
mass m with the resonant fc-mode of an ideal cavity. It is 
2 Interaction of a travelling qubit with a supposed that the atomic flight time inside the cavity is 

, , sufficiently small to treat classically the degree of freedom 

Single cavity mode ^ j a 

of the atomic center of mass along the direction orthogonal 
To pick up the spatial decoherence effect solely caused by to the cavity axis. On the contrary, the atomic transverse 
the entanglement with the qubit-field subsystem, we will dynamics in the x-direction, along the cavity axis, is quan- 
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tized and initially described by a packet of width narrow describes a generic (pure or mixed) state of the cavity 
with respect to the wavelength A of the resonant mode, field, expressed in terms of Fock states The qubit state 
and centered near a nodal region of the sinusoidal mode \ipq{0)) is a coherent superposition 

function. In these conditions it may be approximated by , /„sx 1 \ \ , i6 ■ 1 \ \ n ^ ^ /n\ 

^ IVg(O)) = cos- |e) +e*"^sm- I5) , < 7 < tt (8) 

the linear term, as the interaction part of Eq. ([T]) shows. 

of excited |e) and ground \g) states. 

The conjugate variables x and p just describe this trans- 

We are interested to the behavior of the atomic spatial 

verse dynamics, a and are the usual field operators, and 

coherences, so we consider a coherent superposition of two 

£ is the atom-field coupling constant. Finally, the 1/2 spin 

nearly distinct kets [HHI] , 

operators Sz and S± — Sx i iSy account for the qubit 

dynamics. 1.(0)) = ^[|..(0)) + |..(0))] (9) 



The evolution operator for the Hamiltonian ^ may to describe the initial position distribution of the atomic 

be factorized [2^, for example, in the following form, center of mass along the cavity axis. In particular, we will 

^ ( it ^ assume that the ^-representation of \ifj{0)) {j ~ 1,2) is 

t/(t, 0) = exp <^ — —muN fixX ^ exp <^ - ^-^P^ > 

I J I Zmri J given by Gaussian functions, 

^2 



xexp -flATt^/i^p^ e-' ^, t<To {2) / i \i r ~ xo.j) 



^,(.,0)= exp , (10) 



where To indicates the atomic flight time inside the cavity, ^ j • j i- 1 j • • 

" centered m xq.i and xo,2, respectively, around the origin 

^^'^ of the reference frame which is set in a nodal point of the 

/~~: ehk „ /,^ , , i.T.in^-in^ mode function. The normalization constant of state (1^. 
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We suppose that the initial state of the entire system is 
given by 

p(0) = /5,(0)/5,(0)p/(0) (5) 

where 

pt(0) = 1^(0)) ((^(0)1 , p,{Q) = |v5,(0)) (v5,(0)| (6) 

account for the initial configuration of the atomic external 
and internal degrees of freedom, respectively, while 



1 + exp 



(xo4 - 3^0,2) 
d,Axl 



(11) 



accounts for the eventual (small) overlap of the two Gaus- 
sians. For simplicity, the Gaussian distributions pop are 
of minimum uncertainty, Ax^Ap^ = Ti/2, with the same 
widths Axo and Ap^ for both the Gaussians. 



2.2 Time evolution of the full density operator 

Let us consider the density operator x(0) = /5g(0) ® pf{Q) 
of the irrelevant subsystem. In terms of the dressed states 



pm = Y.^n,n.\n){n'\ (7) |;,±) = i^[|e,n)±|.g,n+l)], |g,0), (12) 
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it assumes the following form 

n,n' —0 

) {Xn'\ +Cn,n' \xt) {Xn'\} 
+ lf2{D^i\^n) + \Xn)){9A 



+En{\xi)-\Xn)){9,0\} + ^\g,0) {g,0\ 



+h.c. 



(13) 



where 



An,n' = ^ {Cn,n' COS^(7/2) + [c„^„' + ie 

+Cn+i,n'e*'^] cos(7/2) sin(7/2) + Cn+i,n'+i sin^ (7/2)}(14) 

Bn,n' = ^ {Cn,n' COS^(7/2) - [c„,„' + ie~'"^ 

+Cn+i,n'e*'^] cos(7/2) sin(7/2) + c„+i^„/+i sin^ (7/2)}(15) 

Cn,n' = {Cn,n' C0S^{l/2) - [c„,„' + ie"*"^ - C„+i,„'e*"^] 

X cos(7/2) sin(7/2) - c„+i,„,+i sin2(7/2)} (16) 

Dn = V2 cos(7/2) sin(7/2)c„,oe~^*, (17) 
En = %/2sin2(7/2) 

Cn+ 1,0, (18) 

F = sin2(7/2)co,o. (19) 



We note that the dressed states (|T2j) are eigenstates of the 
observables N and jlx which appear in the expression Q , 

A.|x±)=±i|x^), N\xt)^{n + l)\xi)- (20) 

Applying the evolution operator ([2]) to the initial state 
/5(0) = pt(0) (g) x(0) and using Eq.s (H]), ((121) and (HOI), we 
get the state of the entire system at time t < Tq, 



X _{<l>n'.lit)\ + {<t>n'.2it)\]<»\Xn){x;A 
+Cn.n' [Kl(i)) + K2W)] 

X _(fc,iW| + ('^,T',2W|] ®|x^)(x,;'|} 

OO 

+ (ICiW) + 1^,2 W)) |X,T)] (.9,01 ((^i(i)l + (¥'2(<)|) 

+i?„ [(Ki(t)) + K2W)) |x,t) " (ICiW) 

+ \K2it)))\Xn)] (.9,0|((^i(i)| + (^2(t)|)} 



4(5 



^ l5,o) (5,01 <E> [\Mt)) + \Mm [{Mt)\ + (Mm 



+h.c. 



where 



(21) 



(22) 



describes the free evolution of the atomic center of mass 
for zero excitations in the qubit-field state, (last term of 
Eq. dni)), while 

it 



(t)) = exp <j ^-mantx | exp <j -^^P' 
X exp<j ±^a„i2pl \ipj{0)) 



(23) 



are the scattered components of the spatial atomic states 
due to the optical SG effect. In fact, a„ ~ ao^/n + 1 ac- 
counts for the acceleration along the cavity axis of these 

components. In the x— representation one has 

1 

Axo 



exp ( T'^™ UnXt 



X exp ■ 



4 Pit) 



(24) 



with 



^(*) = 4J E e-^''''^*''"'"''{^n.»'[Ki(0) <^.(i)=^o,,Ta„iV2, /3(i) = Ax^ + zat/2m. (25) 

n.n'— 

|'/'n,2(^))] [^'?^n',i(^) '^(^'i'n'aW ® \Xn) {Xn'\ For cach valuc of ri wc get au Uniformly accelerated Gaus- 

+Bn.n' [|0n 1 {t}} + |0ri 2(0)] ^^^^ distribution, with the typical width of a free particle. 
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This behavior of the width is a consequence of the hnear- where we have used the fact that the coefficients A„^„, Bn,n 
ity of the positional potential energy which appears in the and F are real. To analyze the atomic spatial coherences, 



interaction term of the Hamiltonian ([T]). 

3 Decoherence Factor 

As said, in our case the irrelevant part is the qubit-field 
subsystem. In other words, we suppose to measure the 
atomic position whatever the values the qubit-field vari- 
ables take on. From a mathematical point of view, this 
implies a tracing of the density operator p(t) on the ficld- 
qubit variables (for a discussion on the relation between 
a reduced density operator and the actual measurement 
of an observable, see [MIH]). In this section we give the 
exact expression of the reduced density matrix that de- 
scribes the atomic translational dynamics. In addition, we 
will see that under sufficiently wide conditions, it is pos- 
sible to single out a decoherence factor. 

3.1 Spatial density operator 

Using the orthonormality of the dressed states (|12p we 
get the following reduced density operator describing the 
atomic translation degree of freedom along the cavity axis. 



Pspace{t) — Trfield,qubit[p{t)] 
-J oo 

-j2{A,,,r.[\cPiM + K,2m 



X [«i(<)|+«2W|]+^« 



we take the matrix elements of Pspace{t) with respect to 
the position eigenstates. Using Eq.s (|22p and and 
separating the real and the imaginary parts of the expo- 
nents we get 

p{x, x'; t) = {x\ Pspace{t) W) 
2 



F 



26V2nAxi{t) 



exp iopQ^ {x^x) 



X exp 



{ 4zix2(t) 



1 



X exp 



1 



^„,„exp iall'^^{x,x') 



+ - 



AAxfit) 
1 



B„,„exp ia^;'l(a;,a;') 



X exp 



25V2nAxi{t)^^^^^ 
1 



AAxfit) 



(27) 



where we have set 



"o [x,x) = 



j-k / i\ 

<,±(a;,a; ) 



8niAxlAxf{t) 

{x - X(i,jf - {x' - XQ^kY 

ht 



(28) 



8mAx^Axf{t) 



■Xt)]' 



1 



^' -<k it) r Tmant{x - x') (29) 



and we have used 



(30) 



\(3{t)Y = Axf^iAx^Q + Apf,t'/m') = Axf^Axfit). (31) 



Once again, we wish to stress that the Eq. (^7)) is the exact 
reduced density operator describing the atomic dynamics 



X [\KM + |0-2(t))] [(0-i(t)| + (C2W|]} 

+ + \'^2{t))] [{•fi{t)\ + {•^2{t)\] , (26) along the cavity axis in the optical SG model, for generic 
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states of the cavity field and the qubit. In the next sec- 
tion the decoherence process hidden in Eq. ([27]) will be p{x,x';t) ~ p{x,x';0)D{x,x';t) (34) 

analyzed for some specific states of the cavity field. , 

where 

p{x,x']0) 



3.2 Decoherence factor ' 2SV2t:Axo 

2 

X 



X1°^P] 4^;^ {x - xojf + {x' - xo^kf [,(35) 



We now observe that the linear approximation of the cav- 
ity mode function requires small displacement and spread 

oo 

of the atomic wave packet, smaU with respect to the wave- ^ ! = 5Z {(^"," + ^n^) cos [k{x - x) ]7„t] 

length A. In other words, when the atom leaves the cavity — i (A„ „ — Bn,n) sin [k {x — x') i7„t]} + F, (36) 



its spatial distribution must be essentially unchanged with 
respect to the initial one. As a consequence, the fiight time 



and 



n + 1 



To must be sufficiently short to allow the disregard of both <">+!' ^ o\/< n > +1. (37) 

the displacement a,,T§/2 and the distortion ApoTo/m in ^ of easy reading and, at the same time, 

Eq.s (El and iB). For the values of the parameters used ^ ^^^^ g^^^ approximation. In fact, all the figures describ- 

in this paper, these conditions are fulfilled for ToH < 10^, .^^ decoherence effects presented below, has been ob- 



where Q - sV< n > +1 is the Rabi frequency. For flight ^^.^^^ j^,^^ ^^.^ equation, but actually they cannot be 
time satisfying this condition we have (see Eq.s m and distinguished from those related to the exact expression 
(EB)) x^,j{t) — xo,j and Axf{t) ~ Ax^. Consequently, 



Eq.s (IMl) and 1^ become Expression ^E^, is suitable also to explicate the relation 

j.k, /x ^ \, \2 j I -.2] /„„x between the field coherence properties and the one-sided 



'^n%{x, x') ~ a{^^{x, x') =F —mant{x — x'). (33) 



atomic defiection [28 in the optical SG effect. A mean 
value of p different from zero requires an imaginary part 
As for a'^^ix, x'), the values of the parameters used in ^^^^^ density operator. For an initial real density operator 
the paper give max {af (a;, x')} < IQ-^ and the phase ^ ^^^^ ^^^-^ -^^^^-^^^ fo^, ^ > 0^ Im[D{x, x'; t)] + 0, 

factors in the second and third terms of Eq. m wiU only ^^^^ ^^^^ ^ ^^^^ ^^^^^ ^^^^ djgjeult to show that 
depend on mant^x — x')/h which takes into account the oo 

<p>{t)^hkJ2 {Bn,n - Au„) ^2nt. (38) 

optical SG effects in the atomic momentum distribution. 

Taking into account all these approximations, the Eq. Looking at Eq.s and we conclude that the atomic 
(j27p may be written as the product of the initial spatial center of mass may acquire an increasing momentum along 
density matrix and a decoherence factor, the cavity axis if its qubit interacts with a field that owns 
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Fig. 1. This figure siiows tlie reduced density matrix p{x, x' \ 0) 
of Eq. (|35p and represents tlie initial condition for all the cases 
considered below. It describes the initial atomic populations 
(near the x = x' direction) and coherences (near the x — —x' 
direction) for the position degree of freedom along the cavity 
axis. The values of the parameters are: m — 10~^^ kg, £ — 
10^sec-\ A = 0.6 X 10"^ meters, Axq = A/100, xqa = -A/20, 
xo.2 = A/20. The positions x, x' are considered in units of jso.il. 

coherence properties, that is, when Cn,n' is not propor- 
tional to Sn,n' ■ 

4 Some examples of field statistics 

In this section the behavior of the atomic spatial coher- 
ences will be studied considering some examples of field 
states. To this end, it is sufficient to specialize the co- 
efficients An,n,Bn.n and F of Eq.s (|14m9p . according to 
the specific initial field state, and then use Eq. ([36]) . The 
spatial initial condition for all the cases here considered is 
given by p{x, x'; 0) of Eq. (j35p . and shown in Fig.[TJ Under 
appropriate conditions [2J, the two peaks near the x — x' 
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direction evolve towards the atomic populations (the two 
possible atomic positions) for the translational degree of 
freedom along the cavity axis, while the peaks near the 
orthogonal direction, x — —a;', undergo a quenching when 
the system evolves towards a classical behavior. 

4.1 Incoherent field states 

From our particular point of view, a common interest- 
ing feature of the incoherent field states is the equality, 
for each n, of the coefficients A„.„ and -B„.„ in Eq. (|36p . 
that implies the reality of D(x,x';t) and, in our case, of 
p{x,x';t). 

4.1.1 Thermal state 

A cavity field in a thermal state at the temperature T is 
described by the normalized Boltzmann factor, 

oo 

= (l - e-'^'''"^) J2 e-"'^^^'^ \n) {n\ (39) 

n=0 

where [3b — ^/ksT and ks is the Boltzmann constant. 
Comparison with Eq. ([7]) gives 

c„,„' - (1 - e-^-''") e-"'5-''"<5„,„. . (40) 

Consequently, for the coefficients A„.„, Bn,n and F of Eq.s 
dUl), (Uni) and HH) we have 

A„,„=iJ„,„4(l-e-^«-)e-"^^- 

X {cos2(7/2) + e-''«'"''^sin2(7/2)} (41) 
F==sin2(7/2)(l-e-'^«'*'^). (42) 
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As said, the equality A„^„ — B„^n implies that the deco- 
herence factor, 

D{x, x';t)= sin2(7/2) (l - g-'^^''") 

+ (1 - e->^^^^) [cos2(7/2) + e-^"''^ sm\j/2)] 

OG 

x^cos[k{x- x')nnt]e-"'^''''"^. (43) 

is a real function. The corresponding density matrix p{x, x'] t) 
= p{x,x';0) D{x,x';t) is shown in Fig. [51 The entangle- 
ment with the field-qubit environment causes a spatial 
coherence decay (upper), followed by a very partial re- 
vival (lower), because of the discreteness of the states to 
which the environment can accede. It is evident that by 
increasing T the revivals should be postponed in time. 

One may also analyze the coherence behavior by look- 
ing at the fimction p{x,x';t) along the x' = —x direc- 
tion. Fig. [3] shows a temporal sequence of the function 
p{x' = — x; t) for the same thermal case. In what follows, 
the spatial coherences will be only analyzed in terms of 

function of this kind. Fig. 2. Position density operator p{x,x';t) of Eq. (|34|) at time 

t = 100 J7~^ (upper) and t = 1000 (lower), when the 

4 12 Coherent state with random phase ^^^'^ thermal state pQ)! . and the decoherence factor 

D{x,x';t) is given by Eq. (|43|l . The parameter values of the 

A mixture of coherent states with random phases is given cavity and the qubit are: T = 200° A", corresponding to < n > 




by 



: 82.76 and 7 = n/2. The other parameters are as in Fig.l. 



P/(0) 



1 

2^ 



„2n 



Z-^ r?,l 



n=0 



(44) 



and using Eq.s ()14m9p we obtain 



In this case we have 



„2n 



(45) 



|2n 



2nl 



„2 7 , |a| 



^1= 



sin' (7/2) 



a| .27 
sm — 

71+1 2 



(46) 
(47) 



10 



A. Vaglica and G. Vetri; Spatial decoherence factor 



The decoherence function, 



Dix, x'; t) = e-l"l' sin2(7/2) + e"!"!' ^ 



7,1 



X I cos^ I + sin^ I ) cos [k {x - x') Q^A (48) 



is equal to real part of the same function relative to the co- 
herent state (see sec. 4.2.1). Consequently, Fig.[S]describes 
also this case. 

4.1.3 Fock state 

Finally, for a Fock state 

Pf{Q) = |no) (nol , (49) 
it is Cn^n' = 5n,no5n' ,no- Consequently, one has 

1 2 7 r -27 

o«,„oCos -+0„,„(,-ism - 



F = 5n,fi sin2(7/2) 



(50) 
(51) 



and 



D{x, x'\ t) = cos^(7/2) cos [k {x - x') Hngt] 
+ sin2(7/2) cos [k {x - x') Hno-it] + sin2(7/2)5„„,o.(52) 

The behavior of p(x' — ~x; t) for the number case, shown 
in Fig. [31 outlines how the coherences survive at all times if 
the exchange of energy (and information) between the sys- 
tem of interest and the environment follows a few paths. 

4.2 Field states with phase properties 

When the field density matrix Cn,n' — Pfi^) of 
Eq. ([7]) is not diagonal, the coefficients A„^„ and i3„,„ 
of Eq. may be different between them. D{x,x';t) and 
p(x,x';t) are in general complex functions. 
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Fig. 3. This figure describes in a simple way the decoherence 
efi^ect (followed by a partial revival, curve (f)), for the thermal 
state (|39p . It represents the function p{x,x';t) defined by Eq.s 
(pi)) . ((35)) and ((43)) . along the x' = -x direction. From (a) 
to (f) time is 0, 3, 10, 100, 200, 500, in units of J7~\ The other 
parameters are as in Fig. 2 

4.2.1 Coherent state 

For an initial coherent state of the cavity field, 

a"(a*)"' , 



P/(0) = e-H Y: rTTJ 
, „ Vn'.n'i 



we have 



and 



V n\n'\ 



a = \a\e 



a 



2n 



2nl 

2n 



7 ae*(''+'^) 
cos ^ + ' ' 



• 7 
sm — 
2 



Sri 



2n\ 



7 |a|e^(«+'^) . 7 

cos — , sm — 

2 V^TTT 2 



(53) 
(54) 

2 

(55) 
, (56) 



while coefficient F is still given by Eq. ((47)) . The decoher- 
ence factor is 

D{x, x'- 1) = e-l"l' sin2(7/2) + e"!"!' ^ ^ 



Tl = 
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X cos + 



2 7 , ^3i,2 7 

n+1 2 



E 



-2i|a|e sin7Cos( 

A2n 



iWn + 1 



sin [k (x — x') fint] ■ 



X') flnt] 




A 




-2 -1 


1 2 


(57) 


-1 





whose real part coincides with the decoherence function 
We recall that, for < n > sufficiently large (as it is 
in our case) and 7 — 7r/2, the phase relation between the 
qubit and the field may lead to two opposite behaviors. 
For 9 + (j) = 7r/2 one has A„.„ ~ Bn^n, D(x,x';t) and 
p{x,x';t) are real, and < p > (t) = 0. On the contrary, 
tor 6 + 4> = 0, or tt the decoherence factor is complex 
and a one-sided deflection characterizes the optical SG 



b) 



c) 







Fig. 4. Same quantity of Fig. 3 when the cavity field is in 

the Fock state (|49|) . with no = 83, and 7 — n/2. From (a) 

effect. In fact one has Bn „ = 0, or An n — 0, respectively ^ ,r\ ^- ■ r 1^ 1^^^ ■ r r^-i rm 

' ■ ' ^ to (f) time IS 0, 5, 10, 50, 100, 200, m units of J7 . The other 



and the atom is scattered to left or to the right (see Eq. 
([551) ). This is the so called trapping condition, after the 
fact that the qubit population is approximately trapped 
to the initial value, with the consequent quenching of the 
Rabi oscillations [5511^ . (For the relation between the 
interference phenomenon of the Rabi oscillations and the 
selective deflection in the optical SG model see Ref. |20]). 
For the case 6 + (j) — and 7 — 7r/2, the behavior of the 
real and imaginary part of p{x' — —x; t) is shown in Fig.s 
[5]and[6l respectively. 



parameters are as in Fig. 1. 

\)^^/^a. In this case the field is 

00 

py(0) = (l-|zp) ^"(^*r'l«)(n'|. (58) 

n,n'— 

and 

c„.„' = (l-|;2nz"(z*)"'. (59) 

The coherent trapping takes on for particular relations 
between the field and the qubit parameters. In fact, it is 
required that 



z = e'^ cot(7/2). 



7r/2 < 7 < TT, (60) 



4.2.2 Phase state 

The exact coherent trapping condition pTU^ may be 
achieved when the cavity field is in an eigenstate of the 
Susskind-Glogower phase operator p51[?7] exp [iff) = {n + 



and 9 + cf) = or TT. 

By using Eq.s (|14j|19p and (H^l) we obtain 
1 



-(l-|z|^)|z|^"{cos^(7/2) 



-2\z\ cos{9 + cj)) cos(7/2) sin(7/2) + \z\^ sin2(7/2)} (61) 



-(l-|z|^)|z|^"{cos^(7/2) 
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Fig. 5. Real part of p{x' = —x-jt) when the cavity field is 
in the coherent state 1)53^ with |ap n >= 82.76, and the 
field-qubit system is in the trapping configuration, 6 + (j) — 
and 7 = 7r/2. From (a) to (f) time is 0, 10, 50, 100, 500, 1200, in 
units of J7~^. The other parameters are as in Fig. 1. This figure 
describes also the behavior of p{x' — ~x; t) for the coherent 
state with random phase, see Eq.s (|48)) and (f57|l . 

-2\z\ cos{e + (f) cos(7/2) sin(7/2) -I- \z\^ sin2(7/2)} (62) 
i^ = sin2(7/2)(l-|zp) (63) 

As one can easily see, for + </> = (0 + </> = tt) one 
has Bn^n — {An^n — 0). The previous consideration on 
the interplay between coherent trapping, quenching of the 
Rabi oscillations and selective atomic deflection in this 
case hold exactly. It is also to note that the photon distri- 
bution of state ([58| follows a geometric law. 



P„ = (n|/5;(0)|n)-(l-|zn|z|^ 



< n >' 



(1+ < n >)"+! 

(64) 



Fig. 6. Imaginary part of p{x' — —x;t) for the same configu- 
ration of Fig. 5. This figure shows the growing (from zero) of 
the imaginary part of the atomic position coherence, followed 
by a loss and a partial revival. From (a) to (f) time, in units 
of J7~^, is 0.3, 3, 10, 50, 100, 1200. The other parameters are as 
in Fig. 5 



with 



< n >= 



cot^ (7/2) 



(65) 



(l-|zP) (l-cot2(7/2))- 

The pure state ((58)) owns the same photon statistics of 
the mixed thermal state (1551) . For this reason, the behav- 
ior of Re[p{x' — ~x;t)] is very similar to the thermal 
case and, for the phase state, we report in Fig. [7] only 

Im[p{x' — —x; t)]. 

Comparing the standard deviations of the geometric dis- 
tribution, An ~< n > for large n, and of the Poissonian 
law. An = \/ < n >, we may conclude that for the same 
value of < 71 > the interaction with a thermal (or phase) 
state takes place following more alternatives with respect 
to the coherent case. As a consequence, coherences fall off 
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Fig. 7. Imaginary part of p{x' = —a;; for the phase state (|58|l 
in the trapping configuration {9 + (j} = 0), with 7 = 0.501 91157r, 
for which it is < n >= 82.76. The other parameters are as in 
Fig. 6 
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Fig. 8. Real part of p{x' = — a;-f 2a;o,i; for the coherent state 
with |ap = 82.76, in the trapping configuration, 6 + (f) = 0) and 
7 = 7r/2. We recaU that the position variable x is considered in 
units of |a;o,i|- From (a) to (f) time is 0, 10, 50, 100, 500, 1000, 
in units of J7^^. The other parameters are as in Fig. 1 
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Fig. 9. Imaginary part of p{x' = ~x + 2xQ^i;t) for the coherent 
state of Fig. 6. From (a) to (f) time is 1, 10, 50, 100, 500, 1000, 
in units of Sl~^. The other parameters are as in Fig. 8 

more quickly in the first case, as one can see comparing 
Fig.s[3]and[5]for the real part of the density matrix, and 
Fig.s[7]and[B]for the imaginary part. 

It is to notice that the coherence revivals of Fig.s [51 
[S] and [7] are displayed for times which are borderline with 
respect to the validity of the model. 

5 Decoherence effects on local coherences 

Let us come back briefly to the link between the imag- 
inary part of the spatial density matrix and the atomic 
momentum. It may be interesting to note that this link 
involves some constraints to the decoherence process of 
the imaginary part of p{x,x';t), at least in some regions 
of of the plane (x, x'). Consider again FiglTJ The two peaks 
along the x' — —x direction describe spatial coherences at 
the initial time (which we may denote as nonlocal). As 
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we have seen, the entanglement with the field-qubit en- ences are not subjected to an asymptotical decay, as Fig. 
vironment causes a decoherence of p{x,x';t) in terms of |9] shows. Since < p > (0) = 0, the function Im[p{x' = 

a time decay of these two peaks. Because of the discrete- —x + 2x()^i;t)] grows up from zero, and for t sufficiently 

ness of the environment states density, the considerations large, its amplitude behaves similarly to the real case of 

brought in the previous section only indicate a propensity Fig. [8l with an essential difference: at the accumulation 

to a classical landing, as the figures of the previous section point a;o,i, Im[p{x' = —x + 2a;o,i; t)] is zero, with an in- 

show. creasing slope which accounts for the atomic acceleration. 



As for the other two peaks located in the orthogonal 
direction, one could say that they describe the atomic pop- 



6 Conclusions 



ulations, that is, the two possible atomic initial position mi i- i n ^ ^ i i • r i ^ i r 

ihe optical Stern-Gerlach model is a useful tool for ana- 

[2]. Actually, when x' is near to but not coincident with xjji ^ix- 

*— ' lyze disentanglement and decoherence processes that m- 

X, the function p(x,x':0) describes coherences (which we i . j_ -^i r i rr j u r 

' r\TT/ \ voive environments with a few degrees of freedom both for 

may denote as local) . An intriguing question concerns the t ^ i i- • i i • ^ i- n r 

' a o ^ discrete and continuous variables. It consists essentially of 

time evolution of these two peaks. First of all, the asymp- t n • ^ i i ^ i • i ^ i -^i 

'- ' ^ a Jaynes-Lummmgs two level atom which entangles with 

totic form of the real part should give the genuine popula- , i i j • i • tt i 

'- ° ° ^ its translationai dynamics along the cavity axis. Under 

tions, however, for continuity reasons, the local coherences , , . ^ i r t ■ ^ 

the environmental action ot the Jaynes-Cummmgs qubit 

cannot abruptly vanish. The Fig. [5] reports i?e[p(x, a;'; i)l 1^111 r ^ • ^ 

^ o 1— 1 f Lf-\ I I yj .^yg analyze the decoherence process of the atomic trans- 

along a direction parallel to (x' — —x), thought one of , ,. , , r ■ , ^ ^ht n j j-t j- ^^ 

° ^ ^ " o lation, considered as system oi interest. We find that the 

the two maxima located in xn , . One may note how xnij, r , , 

decoherence features may be encoded into a decoherence 

asymptotically behaves as an accumulation point for the r t^/ / ,\ t r ^ r ^^ ^ ^ 

J J f function (a;, 2; , i). In fact, for the more general one mode 

local coherences. . , t • j\ -c^i. c u ^ • 1 

state (pure or mixed) of the cavity field, the atomic spatial 

As far as the imaginary part is concerned, also in this density matrix may be factorized as p{x^ x'] t) — p{x, x'] 0) 
case a decay similar to what happens for the nonlocal D{x,x';t), where the decoherence factor D(x,x';t) de- 
coherences cannot be granted. In fact, as seen in sec. 3, pends on the statistics and on the phase properties of 
the mean value of the atomic momentum is connected the field. The field statistics affects Z)(x, a;'; i) through the 
to the Im[p{x, x' ;t)]. More precisely, since in our case density of quantum paths followed in the interchange of 
Im[p{x, x' ;0)] — (and < p > (0) = 0), an increasing energy and information between the subsystems, while the 
mean value of < p > requires a non vanishing imagi- phase properties of the field (more precisely, the qubit- 
nary decoherence factor, as one may see from Eq.s (|38l) field phase relations) are related to the imaginary part of 
and ([55)1 . Wc find, in fact, that the local imaginary coher- D{x, x'; t) and to the atomic kinematics. 
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We find useful to distinguish the local coherences p{x' — 
—X + 2x0, n t) in the neighborhoods of the two population 
peaks, X ~ xqa ~ (i = 1,2), from the non local ones 
p{x' = —x;t), along the direction x' — ~x (see Fig. 1). 
As expected, p{x' — —x; t) asymptotically vanishes, with 
decay times depending on the environment states density, 
and consequently, on the field statistics (because of the 
discreteness of these states, the non local coherences are 
subjected to a partial revival at larger times). On the con- 
trary, the local coherences cannot abruptly vanish, since 
their real part has to ensure the asymptotic survival of 
the populations. Surprisingly, also the imaginary local co- 
herences do survive for increasing values of the atomic 
momentum, and the neighborhoods of the two popula- 
tion peaks asymptotically became a sort of accumulation 
points for the local coherences. In our opinion, these pe- 
culiarities of the local coherences deserve further analysis. 
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